Mesoscopic fluctuations of tunneling through double quantum dots 
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We study fluctuations of conductance of two connected in series dots in the Coulomb blockade 
regime. The pattern of the fluctuations turns out to be extremely sensitive to a magnetic field. These 
conductance fluctuations also provide information about the mesoscopic fluctuations of charge of a 
partially-opened single quantum dot, which are hard to measure directly. 
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Charge of a mesoscopic conductor attached to leads by 
high-resistance junctions is quantized. This quantization 
(Coulomb blockade) is progressively smeared with the 
reduction of the junction resistance to the level ~ h/e 2 . 
Most directly, the smearing can be studied by a mea- 
surement of charge or differential capacitance of a quan- 
tum dot connected to a particle reservoir by a tunable 
junction of conductance Gi. This proved to be a hard 
experiment however, performed only very recently [jjj. 

Spatial quantization of the electron motion in a dot 
leads to mesoscopic fluctuations of its transport ||] 
and thermodynamic characteristics in the Coulomb 
blockade regime. Fluctuations of the conductance in a 
single dot in both peaks and valleys of the Coulomb 
blockade are of the order of the average conductance, 
and were studied experimentally in Refs. Fluctua- 
tions of the ground state energy of a partially opened dot 
result in a random contribution to its charge and differ- 
ential capacitance. However these fluctuations are small 
j|, therefore hard to observe directly. The accuracy of 
the experiments [|J , where the average differential capac- 
itance was studied, is insufficient for quantitative study 
of the fluctuations. 

In this paper we demonstrate that a number of features 
of the mesoscopic fluctuations in the Coulomb blockade 
regime is much easier to observe in a double-dot geome- 
try, shown in Fig. [jj The dots are statistically identical, 
with average level spacing A, and charging energy Ec- 
The temperature T is much smaller than A but exceeds 
the level widths associated with electron escape to the 
leads. Conductances of the three junctions satisfy re- 
lation e 2 /h > Gi ^> Gl, Gr. The dependence of the 
conductance G of the system on V g i and V g % in this case 
reveals sharp twin peaks [§. 

We calculate the variance of the peak splitting and 
show that it gives a direct measure of the random com- 
ponent of the ground state energy and charge of a par- 
tially opened dot. The positions of the peak doublets 
with respect to each other also fluctuate. In the ab- 
sence of the fluctuations, the peaks would form a regu- 
lar pattern (see Fig. ^Ja). Random displacements of peak 
doublets manifest themselves most prominently in strong 
fluctuations of the conductance measured along the line 
Vgi = Vg2 = V g on the plane {V g \,V g 2). These fluctu- 



ations are unusually sensitive to a magnetic field. The 
corresponding correlation field is parametrically smaller 
than the field leading to crossover from orthogonal to 
unitary ensemble. 




FIG. 1. Schematic view of the double-dot system. The 
dots are formed by applying a negative voltage to the gates. 
The conductance Gl and Gr of the channels connecting the 
double-dot to the electron reservoirs is controlled by Vl and 
Vr; voltage Vi controls the conductance Gi of the interdot 
channel. Voltages V g i and V g 2 are used to change the number 
of electrons in the dots. 

Qualitative consideration and the main results — The 
dependence of conductance G of the double-dot system 
on the gate voltages (Vgi, V^) forms a three-dimensional 
landscape. The ridges on this landscape, shown by solid 
lines on the map Fig. ||a, correspond to the activation- 
less hopping of an electron on and off the double dot. 
When Gi — > 0, the ridges form a rectangular pattern 
(dashed lines on Fig. ^a). The conductance has peaks 
when both dots are at resonance (crossings of the dashed 
lines). At these special values of the gate voltages, the 
four states with n or n + 1 electrons in the left dot and m 
or m + 1 electrons in the right dot have the same energy. 
Finite inter-dot tunneling (Gi ^ 0) removes this degen- 
eracy. The two states with the total number of electrons 
n + m + 1 hybridize, lowering the energy of the ground 
state by U, with (U) = (In 2/TT 2 ){hG l /e 2 )E c §. This 
results in anticrossing of the conductance ridges and for- 
mation of twin peak structures shown in Fig. |^. The peak 
splitting (see Fig. ^p) in dimensionless units Af = C g V g /e 
is 
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where C g and C are the capacitance to the gate and the 
total capacitance of a single dot respectively, and (...) 
denotes the ensemble averaging. 
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FIG. 2. (a) The map of the conductance of the double-dot 
system (after Q). The conductance ridges are shown with 
solid lines; the width of a line increases with the ridge height, 
(b) The structure of a peak doublet, displaced from the line 
Mi — N<2 due to the ground state energy fluctuations. 

In random dots, the inter-dot tunneling amplitudes 
fluctuate, giving rise to a random contribution SU to 
the energy U. Let us estimate SU. The tunneling cor- 
rection to the ground state energy results from the vir- 
tual electron transitions between the dots. The number 
of one-electron states contributing to the correction is 
~ Ec/A ||]. The contribution of each state fluctuates, 
with the magnitude of the order of its average, because 
of the fluctuations of the wave functions at the point of 
the inter-dot contact. As the result, 
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However, in a single-dot geometry, this quantity is acces- 
sible only through a difficult measurement of the differ- 
ential capacitance. In the structure of Fig. [j], the meso- 
scopic fluctuations of the ground state energy of a par- 
tially opened quantum dot lead to the fluctuations of 
peak splitting, 
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Here /3 — 1 or 2 for the orthogonal or unitary ensemble 
respectively. The numerical coefficient in Eq. (Q) is ob- 
tained in Quantitative description. Thus the double-dot 
layout allows one to study the fluctuations of the ground 
state energy by means of relatively simple transport mea- 
surements. 

We consider also the random displacements of peak 
doublets with respect to each other. In an ideal struc- 
ture, peak doublets form a square lattice, and the line 
Mi = M2 traverses exactly through the apexes of the split 
conductance peaks. In random dots, however, the ener- 
gies E(n, m) fluctuate, so these peak doublets are ran- 
domly displaced away from the line Mi — Mi (Fig. 
We will show that typically the displacement £ of a peak 



exceeds considerably its width. It leads to a wide distri- 
bution of the heights G max of maxima in the conductance 
measured along this line, G(M) = G(M,M). (Note the 
difference in definitions of G max and the height of a peak 
in G(Mi,M 2 )). 

The displacement £ is a sum of two independent terms. 
One of them, £\ oc SU, stems from the fluctuations in 
the interdot interaction. The second one, £2 ~ A, is due 
to the randomness of one-electron spectra in the dots 
B. The peak width in the direction perpendicular to 
the line Mi — M2 is determined by the hybridization be- 
tween the discrete levels mediating the tunneling through 
the double-dot, and is of the order of y/hGi/e 2 A, which 

The ratio of the 



is much smaller than y/ (£ 2 ) -r \c- 2/ 
two latter parameters yields the probability for the line 
Mi = M2 to "meet" a peak, and thus describes the fluc- 
tuations Of G m ax: 
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The exact result obtained in Quantitative discussion dif- 
fers from the above estimate only by numerical factors: 
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Here A w 0.45 and A w 0.23 for the orthogonal and 
unitary ensembles respectively. We see that the variance 
of G max is parametrically larger than (G max ), in all the 
considered domain Gi % e 2 /h. 

The characteristic magnetic field altering the pattern 
of fluctuations in G(M) also can be obtained from the 
above qualitative consideration. The applied magnetic 
field shifts the two discrete levels mediating the conduc- 
tance [jl0|| . A shift exceeding the peak width leads to 
a considerable suppression of a maximum conductance 
G m ax ■ Such a shift can be produced by a small variation 
of the magnetic field, 
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Here S and Et are the area and the Thouless energy 
M of a single dot, <&o is the flux quantum. Field B c 
is much smaller than the characteristic magnetic field 
B a = (Qq/S) y^A /AttEt changing the electron states 
in a single dot |lf[]. Therefore B c docs not lead to the 
crossover from the orthogonal to unitary ensemble. 

Quantitative description — To calculate the conduc- 
tance of the double-dot system, we use the Hamiltonian 
in the following form 



H — Hd+ He + Hi e ad + H t , 
H d = Ys£k a l a k + J2tt a \ a . 
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The subscripts k and Z, and p and q denote the states in 
the left and right dot, and left and right lead respectively; 

are the one-electron energy spectra of the dots, til,r 
are the numbers of electrons in the left and right dot. 
The dimensionless gate voltages are Mi = C g V g i/e, and 
Ec = e 2 /2C is the charging energy of each dot; C g and C 
are the capacitance to the gate and the total capacitance 
of a single dot. The effects of the interdot capacitance are 
neglected. We also neglect the capacitance of a dot to the 
gate of the other dot, which just skews the map of Fig. ||. 
The Hamiltonian of isolated leads is denoted Hi ea d . The 
conductance of the interdot junction is related to the ma- 
trix elements t kl by G, = [2e 2 /h)(\t kl | 2 )/(2ttA) 2 . These 
matrix elements are proportional to the electron wave 
functions ipk,i in the left and right dot at the point of the 
interdot contact r,: 
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The matrix elements tk P and ti q obey similar relations. 
Using the Hamiltonian (||)-(||) we find the conductance 
G of the system. The randomness of the one-electron 
spectra and wave functions results in the fluctuations of 
the heights and positions of the conductance peaks. 

First we consider the mesoscopic fluctuations of the 
conductance peak splitting M s , which is related to U by 
Eq. (g). To find U we have to calculate the corrections 
to the energy of the four participating states near the 
point of their degeneracy, since U = \E^ 2 '(n + l,m) + 
E^ (n, m + 1)] — [E^ (n, m) + EW(n + l,m+ 1)]. As 
long as the interdot conductance G.; is small, the correc- 
tions can be calculated by the second-order perturbation 
theory in tjy . As an example we give here the correc- 
tion only to E(n, m): 



E^(n, m ) = J2J2 
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The Heaviside function here appears as a combination of 
the Fermi filling factors at zero temperature. The meso- 
scopic fluctuations in U occur mainly due to the fluc- 
tuations of the one-electron wave functions ipk,i(ji), see 
Eq. (bh . The fluctuations of £k,i yield a smaller contribu- 
tion H . Using the Porter-Thomas statistics for the wave 
functions ipk I we arrive at 
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which leads to Eq. (0). The corrections ( [To| ) to the 
ground state energy are given by a sum over the one- 
electron levels with energies < Ec ■ Since there are many 
statistically independent terms in this sum, the distribu- 
tion of the fluctuations SU is Gaussian. 

Now we will consider the fluctuations of conductance 
G(M) = G(M,M). The height of its maxima is a ran- 
dom quantity, depending not only on the one-electron 
wave functions, but also on the random displacement £ 
of a peak from the line M\ = M 2 - If £ < A, the electron 
transport occurs by resonant tunneling via two states, 
fco, io- The height of the conductance maximum is then 
given by 
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with k — kg, I = Iq. Here T k oc ^(fl)! 2 and T/ oc 
|^( r -R) | 2 are the widths of the states k and I with re- 
spect to the electron tunneling into the adjacent lead; 
r^, r# are the points of the dot-lead contacts. At £ > A 
there can be no inter-dot resonance on the line Mi = M 2 , 
and the conduction process inevitably involves cotun- 
neling [H through one of the dots. At low tempera- 
ture T < y/EcA only the elastic cotunneling is possible. 
The resulting conductance G max ~ (A/£)GiGL,R yields 
parametrically small contributions to the ensemble av- 
erages. To calculate (Gmax) and (G^^), one needs to 
use Eq. ( ft2"f ) and average over all the random quantities 
entering in it. The statistical properties of tpf.,1 are well- 
known, and now we discuss those of £. 

One source of fluctuations of £ is the randomness of 
the interdot tunneling matrix elements, which leads to 
random shifts of energies of electron states, see Eqs. (||), 
(|l0|). This contribution yields: 



(£ 2 ) = (SU 2 )/2 
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Another random contribution, £ 2 , comes from the fluc- 
tuations of the positions of the individual levels in the 
dots. Its averaging is not quite trivial. In a real experi- 
ment, the ensemble averaging (. . .) of peak characteristics 
is replaced by averaging over a sequence of conductance 
peaks, and/or over a range of magnetic field B. To study 
fluctuations within the orthogonal ensemble, one must 
avoid using B for acquiring the statistics of fluctuations. 
Averaging over the sequence of maxima along the line 
Mi = M 2 is not actually a well-defined procedure, be- 
cause £ 2 is not a self-averaging quantity. Its average di- 
verges as the number N of conductance peaks taken for 
averaging, increases, ([£ 2 (N)] 2 ) oc IniV ||. In this work 
we refer to the following scheme of averaging: We choose 
an arbitrary conductance doublet with some coordinates 
{Mi°\ M^ ) and then study N consequent conductance 
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maxima along the line (AA 1 (0) +AA,AA 2 (0) +AA). The num- 
ber of maxima N should be large enough to exceed the 
"correlation length" the ground state energy. In a par- 
tially open quantum dot, the characteristic correlation 
length of the fluctuations of the ground state energy at 
different gate voltages is ~ N corr — Ec/A {§■ Thus, one 
must average over N ~ N corr peaks which gives 
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The distribution of £2 is Gaussian, except for the tails 
that are unimportant for our purposes. 

The distribution of the energy mismatch £ is Gaussian; 
its variance is (£ 2 ) = (£f) + (£%), since the two contribu- 
tions are independent. The value of £ is determined by 
a large number of one-electron states. Therefore one can 
neglect the correlations between £ and the one-electron 
wave functions of two states fco and lo, involved in Eq. 
(pi). Averaging of Eq. @ yields: 
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The conductances of the dot-lead contacts for simplicity 
here are taken equal: Gl = Gr = G . The numerical 
factors are ensemble-dependent; A\ = \/^/2tt 2 (GOE) or 
Ai = tt 2 /64v^ (GUE), and 
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Substitution of (£ 2 ) [Eqs. @, @] into Eqs. @, (0) 
gives Eq. (g). 

The fluctuations of G ma x are very sensitive to the mag- 
netic field because the leading contribution to (G^J) 
comes from the peaks, whose narrow resonant region is 
traversed by the line A/i = A/2 . A variation of the mag- 
netic field by B c [Eq. (Q)] would shift |l(J these peaks 
away from the line by ~ y / hGi/e 2 A. This leads to a sig- 
nificant decrease of the height of the corresponding max- 
imum in G(Af). The one-electron wave functions are not 
changed by such a small variation of the magnetic field. 
Therefore the correlation function of the peak conduc- 
tance at different magnetic fields, (G pea k(Bi)G max (B2)) , 
is controlled by the parametric correlations of £ 1 with 
itself at two diffe rent fields. If \B\ - B 2 \ < B A = 
($o/<S')\/A/47r£ , T, the distribution function of the field- 
induced variation e(Bi — B2) of the energy £ is 10 
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Using Eqs 
tion (G 
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(p^|)-(p^), (|T7|), we obtain the correlation func- 
:{Bi)G ma x{B2)) ■ Its full-width-half-maximum 
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The full expression for (G max (Bi)G rnax {B2)) is too cum- 
bersome to present here. 

Conclusion — In this paper mesoscopic fluctuations 
of Coulomb blockade in a double-dot structure are stud- 
ied. Electron tunneling between the dots leads to the 
formation of peak doublets in the dependence of conduc- 
tance on gate voltages V g \, V S 2- We find the variance of 
spacings within the doublets, see Eq. (||). This variance 
provides information on the mesoscopic fluctuations of 
the ground state energy of a partially opened quantum 
dot, which is hard to measure directly. The pattern of 
the conductance as a function of V g i, V g 2 is affected by 
an unusually small magnetic field given by Eq. (|J). 
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